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Abstract 

TPEcalc is a C++ program for calculation of two-photon exchange (TPE) amplitudes in elastic electron- 
hadron scattering, based on the dispersion method. It is a command-line tool which accepts kinematical 
parameters (Q^ and e) as input and returns TPE amplitudes. It can do calculations for proton, neutron and 
pion targets. Any spin-1/2^ or spin-0~ target is supported, if the user supplies all necessary form factor 
parameterizations. This paper describes how to use TPEcalc and outlines underlying theory. The program 
source code can be downloaded from http : //tpe . bitp . kiev . ua/ , 



CD . 1 Introduction 

During last decade, two-photon exchange (TPE) has drawn much attention of both theorists and experimental- 
ists. It was shown that TPE corrections remove the discrepancy between Rosenbluth and polarization methods 
in proton form factors measurements. The TPE corrections have been applied to low-Q^ [T] and high-Q^ [3J [5] 
elastic ep scattering. Many experimental searches for TPE effects are proposed (see e.g. [3] and references 
wherein) and some are already finished [5] . 

The theoretical evaluation of TPE effects was done by several authors [2l|6l[7|. However, at present there is 
no publicly available computer program for the calculation of TPE corrections, though the need for it clearly 
exists. Herewith we present the program TPEcalc, which is capable of calculating TPE amplitudes for the 
elastic ep, en and btt scattering, with a variety of intermediate states (ISs). The calculation is based on the 
dispersion method |7j, see also [B]. It closely follows the procedure, described in Sec. HI of Ref. ^, thus there 
is no detailed discussion here. 

The program is written in C-l — h. It is a command-line tool which accepts kinematical parameters {Q^ 
and s) as input and returns TPE amplitudes. It has such options as customizing output format, calculation 
precision and so on. Moreover, the user can easily change particle masses and form factors or append new ISs, 
without changing the code, since they are recorded in the separate text file, which is read and interpreted by 
the program. 



2 Theoretical background 
2.1 Kinematics 

First, we summarize some relations, pertaining to the kinematics of the elastic electron-hadron scattering. 
Let k (k') be initial (final) electron momentum, and p {p') be the respective momenta of target hadron. The 
transferred momentum is g = p' ~P, the Mandelstam variables are s = (p + fc)^, u = [p — k'Y, and t = = —Q^. 
The electron mass m is usually neglected, and the hadron mass is M . 

The kinematics of the elastic scattering of two particles can be fully described by two parameters. We adopt 
for this purpose and e, the virtual photon polarization parameter: 

e=[l + 2{l + qyAM^) tg\e/2)] , (1) 
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where 6 is laboratory scattering angle. The parameter e can also be expressed as 
where v = s ~ u = {p + p'){k + k'). 

2.2 Definitions of TPE ampHtudes and observables 

Target spin 0~ 

In general, the elastic scattering amplitude contains one form factor (FF) and has the following form 

M = -^u'^^^u{p + p')^F{q^e) (3) 

where a is fine structure constant and u {u') are initial (final) electron spinors. The one photon exchange (OPE) 
amplitude has the same form, except that the form factor F depends on only. Thus we may write 

F{q^e) = F{q^) + 6F{q\e) + 0{a^) (4) 

where 5F is TPE amplitude. The quantity output by the program is "normalized amplitude" 5F/F. 
The TPE correction to the cross-section is obviously 

- = 2Re- (5) 
Target spin 1/2+ 

The elastic scattering amplitude has the following structure 



r 



U (6) 



(where U and U' are initial and final spinors for the target), whereas in OPE approximation F3 is absent, and 

J^i, F2 depend on (f' only. 

Instead of F\, F2, F3, it is convenient to introduce 

Qe = Fi- TF2 + vF2./AM^ 

Gm = Fi+F2+eiyF3/4M^ (7) 
03 = vFs/AM^ 

In OPE approximation Qe and Qm become usual electric and magnetic FFs, thus 

QE{q\e) = GE{q^) + 5gE{q\e) 

GMiq^e) = GM{q^) + 5QM{q\e) (8) 
Qz{q\e) = Sg^iq^s) 

The cross-section is more simply expressed through this set of amplitudes, see below. It is common to normalize 
TPE amplitudes for the electron-proton scattering by the proton magnetic FF, thus the program output is 

SGe/Gm, SGm/Gm, SGs/Gm- 

The correction to unpolarized cross-section will be 



a eR^ + T \ Ge Gm , 

where r = Q^/4M^ and R = Ge/Gm- The following relations may prove useful for polarization transfer 
experiments: the corrections to measured FF ratio Rexp, 

SRexp ^ p^(^ / SQm s{1 - e) 5Qz \ ^^^^ 

Rexp I Ge Gm 1 + s Gm J 

and to longitudinal component of final proton polarization, P; , 
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2.3 Method outline 



In this section wc briefly describe the calculation procedure and the formulae used. 
First, the imaginary part of the scattering amplitude can be written 

^2 



R spin 



(12) 



where p" = p + k — k" , ti — (k" — k)"^, t2 = {k" — fc)^, and AA^f^^ are amplitudes for the excitation of the hadronic 
state R by virtual photon (i.e. T7* — i?, where T stands for the target). They contain on-shell transition FFs 
and are explicitly written down in Appendix VK\ 

The full TPE amplitude is then reconstructed using analyticity property. Finally, the TPE amplitude, say, 
5F, is expressed as 



JF(.)=;^[<5^^(f(.)-<5^^(j^)(- 



(13) 



where 



[l/2 +t(4M2 -t)]2 



E / d'k' 

■1,3 



X 



(fc"2 - m2)(p"2 - Af2) 



tlt2 

fc"2-m2 + p"2-M2 + ^M. ' ^2 , ^, t) 



(14) 



where i^/^^ are FFs for the electromagnetic excitation of the IS R. The quantities a\^'' , A^^^- , Ap^] and ^i^]- 
are polynomials in their arguments and actually depend on the IS R only via its spin-parity. 
Parameterizing transition FFs as a sum of poles 



t — mi 



(15) 



as assuming they decrease at t — ?> cxd fast enough, the above formula may be rewritten as 

d^k" 



SK 



(R) 



box 



tjab 



{ti - ml){t2 - ml) 



(16) 



4f (m2,m2,^,0 A^-;(m2,m2,^,t) A';^>^{mlmi,v,t) 



\iR) 



+ A^^Ami,rni,i.,t) 



(fc"2 -to2)(p"2 -M2) 



fc"2 



p"2 - M2 



3 Program usage 

3.1 General 

The program is invoked as follows 

TPEcalc [options] parfile 

All available options are described in detail in Sec. 13.21 Mandatory parameter parfile specifies the file, 
which contains the information about the target and ISs, such as masses, spin-parities, transition FFs. Three 
predefined parfiles come with the program: proton. par, neutron. par and pion.par, for proton, neutron and 
pion targets respectively (see description in Table [IJ . 

At first, the program reads and verifies particle descriptions from the parfile. If everything is fine, it prints 
a short message, like this 

TPEcalc version 1.00 

Current options: -i= -0= -h= -f=a -n=1000 -c=0.000 -m=0 proton. par 

Loaded 2 particle (s): 

(x) proton: M=0. 938270, JP=l/2+ 

( ) Delta: M=l. 232000, JP=3/2+ 
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A mark on the left of the IS label indicates that the contribution of this IS will be calculated and added to the 
result. The particle masses here (as well as other dimensional parameters like below) are in units of GeV 
(GeV^). 

After that, the program reads and e values, either from the stdin or from the file (if -i option is present) 
and calculates corresponding TPE amplitudes. The amplitudes are output in the table form either to stdout 
or to the file (-o option). The table columns are, 

• for target spin 0~: 

Re^ Im^ 

• for target spin 1/2+: 

Re|2^ Im|2E Re Im Re |2a. Im |Sa- 

tJM (JA7 (JA/ (jA/ 

The definitions of the amplitudes and 8Qe, SQm, ^Gs are given in Sec. 12.21 The above-described output 
format may be overridden with the -f option. 

The elastic contribution (i.e. that of the IS coinciding with the target) is well-known to be infra-red divergent. 
If the elastic contribution is included, standard Mo&Tsai correction ([TT], see also [B]) will be automatically 
subtracted to cancel the divergence. 

3.2 Options 

-i=FILE, input file 

If this option is present, program takes and e values from the file FILE, otherwise from stdin. 



-o=FILE, output file 

If this option is present, resulting TPE amplitudes are written to the file FILE, otherwise to stdout. 



-h=ISl ISn, intermediate states to include 

This option allows to specify which ISs (of those contained in the parfile) will be included in the calculation. 
One should supply comma-separated list of IS labels with no spaces between them. If some IS is not found in 
the parfile, it is ignored and a warning message is printed. In absence of this option, ISs are included or skipped 
according to the defaults in the parfile (see Sec. 13. 3p . 



-f=[qeaf], output format 

By default, the program prints only calculated TPE amplitudes (real and imaginary parts). The -f option 
allows to customize the output. One should supply a combination of letters q,e,a,f , which specify the columns 
of the output table as follows: q — Q^, e = e, a = TPE amplitudes, f ~ elastic FFs of the target. E.g. with 
-f =qea the program will print 

Q2 £ Re^ Im^ Re#ii Im 4^ Re im^S^ 

(jAf LtM (-TA/ Ltm Um (jA/ 



parfile 


IS label 


particle 


FFs, Ref. 


proton. par 


proton 
Delta 


P 

A(1232) 


9^ 


neutron. par 


neutron 


n 


[W 


pion.par 


pion 
rho 
bl 


P 

5i(1235) 


8 
8 
8 



Table 1: 
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with -f =qf no TPE will be calculated and 

will be printed (assuming 1/2+ target). 
-n=NUMBER, calculation precision 

The loop integrals, needed for the calculation, are evaluated semianalytically (using numerical integration as a 
final step). The precision of the latter can be controlled by the -n option. The NUMBER must be positive integer; 
the greater the number, the more precise (but also more time-consuming) the calculation will be. The default 
is -n=1000. 

-c=NUMBER, form factor cut-off 

For the program to work correctly, transition FFs should decrease fast enough as oo. To achieve this, 

one may use -c option. The FFs specified in the parfile will be damped according to 

FiQ') ^ F{Q')^J^ (17) 

where A (in units of GeV) is the number, specified after -c. 

-ni=NUMBER, electron mass 

This option is experimental. It sets non-zero electron mass (in units of GeV). In the current version, the 
only place, where the electron mass is present, is propagator denominators. In any case it is neglected in the 
numerators. 

3.3 Parfile format 

The parfile contains the information about the target and intermediate states, involved in the calculation. It is 
a plain text file, and the user may edit it or create new parfiles adhering to the format, described below. 
A line, which begins with a % sign, is a comment (ignored by the program). 

The very first line in the file contains version number, which must match the program version. This is 
necessary because the parfile format may change in the future versions. Please do not modify this line. 

The rest of the file consists of one or more particle descriptions, according to the following format. The first 
particle in the file is the target, the others are inelastic ISs. Here is an example: 

pion : 1 

M=0. 13957 JP=0- 

7„ this line is a comment 

1. 

0.7194 -1. 

rho : 

M=0.7758 JP=1- 

0.725333333 

0.78259 -0.725333333 

Each particle description begins with a label, followed by a colon. The label is used solely to identify the particle, 
and may be any word, convenient to user. It should consist of alphanumeric characters and must not contain 
spaces or commas. The maximum label length is 15 characters. A digit following the colon is the "inclusion 
flag". If it is non-zero, then the contribution of this IS is calculated and added to the resulting TPE amplitude. 
Otherwise, the IS is skipped. This behaviour may be overridden with the -h option. 

Next line specifies mass (in units of GeV) and spin-parity of the particle. The mass can be any positive 
number, supported spin-parities can be found in Appendix [X\ 



5 



The subsequent lines contain a matrix of the transition FFs (for inelastic IS) or elastic FFs (for elastic IS), 
terminated by a blank line or end-of-file. The matrix has the form 

mi Cll ■•■ Cfcl 

m2 Cl2 ... Cfc2 ^-j^g^ 

and defines the transition FFs via 



Thus the number of columns in the matrix is (number of transition FFs) + 1. The first row should always 
contain zero mass, toi = 0, and the coefficients Cn = Fi(0). The corresponding transition amplitudes in terms 
of FFs Fi for each supported spin- parity are written down in Appendix 

3.4 Limitations 

The input parameters must satisfy > 0, < e < 1. In the low-e region, the accuracy drops seriously as 
e — ^ 0, because the calculation method involves a subtraction at e = 0. Setting higher precision with -n option 
may help. Though e = is legal kinematics (strictly backward scattering), the calculation at e = is impossible 
due to the above property. 

With the inelastic ISs, the accuracy also drops when the total cm. energy squared s is near the resonance, 
s M|. For m = TPE amplitudes diverge as In |s — near that point. This is an artifact, arising from 
neglecting IS width. Actually, if one takes finite width into account, the s dependence become "smeared" and 
the divergence should disappear. However, this is not possible in the present version of the program. 

4 Compiling the program 

We have successfully compiled TPEcalc with GNU gcc (g++), under both Linux and Windows, but compiling 
with any other C-I--I- compiler should constitute no problem. First, unpack TPEcalc . zip into an empty directory. 
You should have 15 source files (of which 5 in the subdirectory A/), 3 parfiles (proton. par, neutron. par, 
pion.par), and makefile. Then run make to build the TPEcalc executable. If you use some other compiler 
(not gcc), change two first lines in the makefile: 

cxx = <name of your C++ compiler> 
options = <compiler options> 

It is good idea to set maximum optimization, to ensure best performance. 

The parfiles should be in the current directory when running TPEcalc, otherwise you must provide full path 
to them. E.g. under Linux you might want to place TPEcalc and parfiles into $HOME/bin, and then invoke the 
program like this 

TPEcalc [options] ~/bin/proton.par 

A Structure of the transition ampUtudes 

This section contains the formulae for electromagnetic transition amplitudes between supported targets and ISs 
(T + 7* R, where T is target and R is IS). Below, p and q are the momenta of the resulting hadron and 
virtual photon, respectively, thus p — q will be the target momentum. Mn is the mass of the IS, M is target 
mass. In each case, the program assumes that the parfile contains corresponding form factors, Fi, F2, etc. 




(19) 
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A.l Target spin 

• For the elastic IS 



- {2p,, - q^) (20) 



• For pseudovector IS (J^ = 1 ), e.g. p meson: 



^^=MTZ7J^^'"'^^1-^P'^^rP^ (21) 



2Mr 

Here and in the next equation v is spin-1 wavefunction of the IS, satisfying Pav" ~ and Vav" = —1. 

• For vector IS {J^ = 1+), e.g. bi meson: 

= j^p^^p (g^^g*^ - a'^'qniPaFi - q^F2)v. (22) 

R 

A. 2 Target spin 1/2+ 

• For the elastic IS 

Mf. = U' (^^^F,^[j^,q]^^U (23) 

i.e. Fi and F2 are usual Dirac and Pauli FFs. Note that TPE amplitudes will be normalized by Gm = 
Fi + F2. 

• For the IS spin-parity 3/2+ (e.g. Delta resonance) 

= ^i9^"l'^ - 9''"<f ) U [{Plu - P.)Fi +p,F2 + quF^] 75 y« (24) 

where is Rarita-Schwinger wavefunction, normalized according to VaV" = — 2Mj^. Here Fi corresponds 
to magnetic, F2 — to electric and F3 — to Coulomb form factor. 
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